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ABSTRACT 

Short period planets are subject to intense energetic irradiations from their stars. It has been shown that this can lead to significant 
atmospheric mass-loss and create smaller mass planets. Here, we analyse whether the evaporation mechanism can affect the orbit of 
planets. The orbital evolution of a planet undergoing evaporation is derived analytically in a very general way. Analytical results are 
then compared with the period distribution of two classes of inner exoplanets: Jupiter-mass planets and Neptune-mass planets. These 
two populations have a very distinct period distribution, with a probability lower than 10~ 4 that they were derived from the same 
parent distribution. We show that mass ejection can generate significant migration with an increase of orbital period that matches 
very well the difference of distribution of the two populations. This would happen if the evaporation emanates from above the hottest 
region of planet surface. Thus, migration induced by evaporation is an important mechanism that cannot be neglected. 
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1. Introduction 

Different scenarios have been proposed to explain the forma- 
tion of hot Nept unes and Super-Earth s: migration in a pro- 
toplanetary disk ( Mordas ini et al.ll2009l). in situ form ation by 

| accretion of planetesimals dBrunini & Cioncol I2005I) , embryo 
formation in a compact multip lanetary system and subsequent 
migration through scattering ( Ida & Lin N2010l) . tidal downsiz- 
ing of giant planet embryo (lNayakshinll2010l l201ll). or par- 
tial evaporation of more massive planets dBaraffe et al.l [2004; 
Lecavelier Des Etangs 2007; Valencia et al. 2010). Here, we fo- 
cus on the evaporation hypothesis with a particular attention to 
semi-major axis evolution. Indeed, among the shortest period 

\ exoplanets, those cl oser to their star tend to be more massive 
than the outer ones (iMazeh et al.ll2005l: [S outh worth e t~ai1l2007l: 
iDavis & Wheatlevl2009l: lBemtez-Llamb av et al.ll201ll) . This has 
been attributed either to a combined eff ect of tidal interactions 
and m igration in the protoplanetary disk ( Ben itez-Llambay et al.l 
l201ll) . or to the fact that low- mass close-in planets can not sur- 
vive catastrophic evaporation (IDavis & W heatley 2009). 

While isotropic planetary evaporation generates almost ne- 
glectable migration, anisotropic ejection of matter can lead to 
substantial increase of semi-major a xis. This effect is simi- 
lar to anisotropic thermal radiation ( Fab ry ckyl 120081) . In the 
case of hot gazeous planets, two different mechanisms have 
been proposed to explain atmospheric escapes, i.e. how par- 
ticules reach the Roche lobe (Hill) radius of their planet. On 
the one hand, this may be due to a thermal inflation of the 
planet. Particules of the upper atmosphere are assumed to fol- 
low a Maxwellian velocity distribu ti on, and those with the high- 
est speed escape dGu et al.l 120031; iLecavelier des Etangs et al.l 



120041) . In that case, the evaporation is preferentially directed to- 
wards t he LI Lagran gian point located between the star and the 
planet dGu et al .1120031) . On the other hand, the es cape may be 
driven by vertical winds on the planet's day side (lYelld [2004; 
iGarcia Mun"ozll2007l: iMurrav-Clav et al.ll2009h . This situation is 
more complex and 3D hydrodynamical simulations are required 
to get the geometry of ejection. Nevertheless, winds are likely 
anisotropic. 

Atmospheric circulation models predict the existence of a jet 
that displaces the hottest regi on of hot Jupiters from the s ubstel- 
lar point, along the equator ( Showman &Guilloti r2002). More 
recent models agree that the jet is a common feature of circula- 
tion models, and that it is confined to latitudes between -30° and 
30° with a displacement from the substellar point usually east- 
ward, by 10-60° of longitude. Its value depends on the strength 
of the imposed stellar heating and other f actors such as the ra- 
diativ e timescale and drag time constant (IShowman & Polvanil 
l201ll) . These predictions are corroborated by the observed sur- 
face temperature distribution of HD 189733b ( Knu tson et al.l 
120071) . 

Although there are evidences of a displacement of the hottest 
region with respect to the substellar point on planet surface, the 
dynamics of exobase is poorly known. Strictly speaking, it is 
thus not possible to assert in which direction and how mass is 
ejected by evaporation. Hereafter, in order to get estimates of the 
migration induced by evaporation, we will take the results men- 
tionned above concerning planetary surface as representative of 
the exobase structure. 

This letter is organized as follows. In Section [2l we derive 
analytically the orbital evolution of a planet undergoing evapo- 
ration in a very general way. In Section O we compare the dis- 
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planet 



axis. Another scaling would modify only slightly the eccentric- 
ity evolution. The mass loss-rate is thus inversely proportional 
to the square of the distance 



m(t) = -B(t)/r 2 



(3) 



where B(t) is any slowly varying function accounting for the 
time evolution of the stellar luminosity. From ©, 0, and the 
expression of V, one gets 



Fig. 1. Geometry of evaporation. 



— = -fi(t)— + K(f) sin0— 



(4) 



tribution of two classes of close-in planets: Jupiter-like planets 
and Neptune-like planets. The statistical analysis is then used, 
in Section ffl to constrain the geometry of evaporation assum- 
ing that hot Neptunes are partially evaporated hot Jupiters. This 
geometry is also confronted with the predictions of atmospheric 
models. The conclusion is presented in Section 

2. Orbital evolution with evaporation 

Two-body problems in which one of the component loses mass 
should be solved with caution. Severa l misanderstandings hav e 
already been found in the litterature dPlastino & Muzzioll 1992b . 
since the equations of motion are not the same whether the ejec- 
tion is isotropic or not. Here, we consider a simple model of 
evaporation that is well suited for both isotropic and anisotropic 
mass-loss. 

Let a planet of mass m(t) orbiting a star with mass m* . The 
orbit is characterized by a radius vector r, an orbital velocity 
v = r, a semi-major axis a, an eccentricity e, an orbital angular 
momentum t - r x v, and a period P = 2n/n, where n is the 
mean motion defined by G(m + m+) = n 2 a 3 . G is the gravita- 
tional constant. Let v p and v* be respectively the velocity of the 
planet and of the star in a fixed reference frame. The orbital ve- 
locity is then v = v p - v*. For any vector w, we note w = \\w\\ 
its norm, and w = w/w its unit vector. The frame (r, 0, k), where 
k = € and 6 = k x r, is used to defined the geometry of the 
mass ejection, which is given by two parameters: the angle be- 
tween the star and the direction of mass ejection, and the aper- 
ture (p of the stream (see Fig.Q]). In particular, the evaporation is 
completely focused in one direction when (p = 0°, and becomes 
isotropic when <p = 180°. We note u = -cos Of - sin 06 the 
average direction of the stream. Let V e s be the ejection speed 
of atmospheric particules with respect to the planet barycenter. 
Note that Ves can vary between 10 and 20 km.s -1 for a planet 
orbiting a Main- Sequenc e star, and a T Tauri star, respectively 
(Murrav-Clav et al. 2009). The mean velocity V of the stream is 
T]((f)V es u, where 



cos <p' sin <p' dip' 
sin <p' dip' 



2 <P 

= cos - 
2 



(1) 



(2) 



and the equations of motion are 

dv p r 
m = -Gm+m— + mV , 

dt r 3 

dv+ ^ r 
rn±—— = +Gm*m— . 

dt r 3 

We then assume that the mass-loss rate, m(t), scales linearly with 
the energy received by the planet from the stellar XEUV emis- 
sion. It should be stressed that this has no effect on semi-major 



where K(t) = rj((p)B(t)V e Jm(t) and ji{t) = G(m(t) + m*) - 
K(t) cos 0. Notice that even in the extreme case with a mass-loss 
rate m = -10 15 g.s -1 , and an ejection speed V es = 100 km.s -1 , 
for a hot Jupiter at 0.05 AU, K/(Gm) remains lower than 2 x 10" 8 . 
In the following, we thus use simply ji(t) = G(m(t) + m*). 
From the definition of the semi-major axis 



2 r 



2a 



and the expression of the Laplace-Runge-Lenz vector e 



e = r . 



(5) 



(6) 



the equation of motion © leads to 



dt \ r/fi 



a 1 Ik sin 
\ir 3 



(7) 



de _ s u KsmO , n 9 * 

— = - (r + e) - + — (fir + rv x k) . 

dt fi jir 5 

The first terms of the right-hand sides of Eq. (|7]) are only due 
to the mass decrease. They do not depend on the geometry and 
remain unchanged even for isotropic evaporation. On the other 
hand, the last terms of Eq. (|7]) vanish for isotropic (k = 0) or 
radial (sin = 0) mass ejection. The long term evolution of the 
orbital parameters is obtained by averaging over the mean 
anomaly M keeping u ,(t), K(t) and fi = na 2 Vl - e 2 constant (e.g. 
iBoue & La skar 2009). The results are 



da \ ii naK sin 

— ) = -a- + 2 



dt \u 
de 
dt 



li ji(\ - e 2 ) ' 
HK sin 

-e . 



(8) 



fi(l + Vl - e 2 ) 



In the following, we consider the simplest case where the orbit is 
circular. Eq. ® shows that such orbits remain circular over the 
whole evolution. Replacing B(t) (see Eq. ©) by -a 2 m(t) in the 
expression of K(t), and using v = na and n(t) = G(m(t) + m*), 
the secular evolution of the semi-major axis becomes 



/// 



M,e=0 



m+ + m 



v m 
■ 2t a . 



(9) 



where vo = v(t = 0) is the initial orbital velocity, r is a dimen- 
sionless parameter that controls the efficiency of the migration. 
We will thus refer to this parameter as the 'migration efficiency' . 
Its expression is 



2 <P ■ n V ™ 

t = cos — sin — 

2 v 



(10) 
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For isotropic evaporation (cp - 180°) or radial ejection (0 = 
0°), the migration efficiency vanishes. In that case, the prod- 
uct a{m+ + m) of the sem i-major axis with the t otal mass of the 
system remains constant ( Hadjidem etrioul [19631) . Thus, even if 
a Jupiter-mass planet orbiting a Sun-like star loses all its mass, 
the semi-major axis increases only by a factor Mj up /M = 0.1%. 
The term amjim* + m) in Eq. © will thus be neglected in the 
following. 

For arbitrary migration efficiency, the solution of Eq. @ is 
a(t) = ao [I + rln(m(t)/mo)]~ 2 . For convenience, we provide 
an equivalent relation, with orbital period instead of semi-major 
axis, deduced from Kepler's third law 



Histogram of Period for all RV planets 



— = (1 + ring) 
A) 



(11) 



where q = m/mo is the fraction of remaining mass in the planet. 
The expression (TTTT) is independent of time, and holds for any 
mass-loss rate evolutions. The only requirement concerns the ef- 
ficiency of evaporation. Here, we assume that the evaporation 
process is sufficiently strong to decrease the mass of a planet 
from mo to m. 

Using a typical ejection speed V es = 15 km.s -1 and an ini- 
tial orbital velocity vo = 145 km.s -1 corresponding to a 3 days 
orbit, the order of magnitude of the migration efficiency given 
by (TTT)b is r ~ Ves/vo = 0.1. With this value, if a Jupiter mass 
planet is transformed into a Neptune mass planet by evaporation 
(q = 0.05), the increase of its orbital period deduced from (fTTb 
is P/Po = 1.5. For comparison, the ratio of the median periods 
of Jupiter-like planets and Neptune-like planets (see section [3]) 
is 4.9/3.4=1.4. 



3. Observed period distributions 

From the more than 700 exoplanets discovered up to date, 204 
have orbital periods smaller than 10 days and mass smaller than 
5 Mj UF Q. One can group these planets in two classes: Neptune- 
mass planets and smaller (m < 2 M^ Qp or 2 x 0.054 Mj up ), and 
Jupiter-mass planets (m > Ms a t or 0.3 Mj up ). The first group 
contains 37 planets, from which we removed two couples be- 
longing to the same system^, while the second 151. Both corre- 
spond to 92% of short-period planets considered. The two popu- 
lations are characterized by very distinct period distributions: the 
Neptunes have a median period of 4.9 days, an average period 
of 5.3 days and a standard deviation of 2.6 days, while for the 
Jupiters the same quantities are of 3.4, 3.5 and 1.6 days, respec- 
tively. A Kolmogorov-Smirnov test shows that the cumulative 
distribution functions (CDF) of the two distributions are signifi- 
cantly different, with the probability that they were derived from 
the same parent distribution being of 3.8xl0" 5 (see Fig.©. 

It is important to keep in mind that the interpretation of these 
results is conditioned by different factors, which warrant a dis- 
cussion. Firstly, the fitted parameters depend on the limits de- 
fined for the planetary classes. However, the general properties 
of the distributions are kept for a wide range of parameters, and 
the results are expected to be resistant to changes of these param- 
eters. More important is that the observational bias imprinted on 
the sample by radial-velocity detections, which depends on mass 



1 Information extracted from Exoplanet Encyclopaedia, 
|http://exoplanet.eu/| on the 29/1 1/201 1 . 

2 For the purpose of this study, we discard all systems containing 
more than one planet with period shorter than 10 days. The removed 
planets are (CoRoT-7b, CoRoT-7c), (G1581b, G1581e), (HD 40307b, 
HD 40307c), and (Kepler- 18b, Kepler- 18c). 




CDF for the two classes of planets 




Fig. 2. (a) Distribution of orbital periods for all planets with 
m < 5 Mj up in white, Jupiters in gray, and Neptunes in black, (b) 
CDFs with their best fit by a log-normal distribution. Jupiters are 
represented by triangles, and Neptunes by circles. 



and period, is different for the two classes of planets. Since the 
amplitude of a signal and thus its detectability is proportional to 
m p-i/3 ? heavier and closer planets are detected more easily than 
lighter far away planets. For this reason, the difference in median 
or average period values is expected to increase for the real unbi- 
ased population, showing that the difference in values we found 
here is not a detection effect but rather an underestimation of the 
real one. In any case, the reflex motion induced by a Neptune on 
a 10-day orbit on a 1 Solar-mass star is of 5 m.s -1 , significantly 
larger tha n the precision of sta te-of-the-art spectrographs (e.g. 
HARPS) (Mayor & Udrv 2008). As a consequence, we expect 
the detections reported up to date to be representative of the true 
population. 



4. Constraints and comparison with atmospheric 
models 

Here, we derive the distribution of the migration efficiency r 
assuming that the hot Neptune and the hot Jupiter populations 
correspond to the same planets at different degrees of evapora- 
tion. This provides constraints on the geometry of evaporation 
which are then compared with results of atmospheric models. 
That analysis relies on the possibility that hot Jupiters can lose 
95% of their mass by evaporation. Although additional works 
are required, we expect this to happen at the early age of p lanet's 
history, when stellar emiss ion is strong dRibas et al.ll2005l) . and 
planet density is low (e.g. lMordasini et al.ll2010h. Indeed, thes e 
two conditions lead to important mass-loss (lErkaev et al.ll2007h . 
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Fig. [3tb)). This value matches very well the one predicted by 
atmospheric models (IShowman & G uillot 20021; iKoskinen et all 
2007 1: Ishowman & Polvanil l201ll) or reported on HD 189733b 
dKnutson et alJl2007h (filled circle in Fig. Eft))). 

Nevertheless, it should be stressed again that the dynami- 
cal models and the observation of the temperature map of HD 
189733b do not extend up to exobase. It is thus not possible to 
prove the connection between the direction of mass ejection and 
the position of the hottest region. 



-0.2 -0.15 -0.1 -0.05 0.05 0.1 0.15 0.2 
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Fig. 3. PDF of the migration efficiency r as a function of r (a), 
or as a function of the angles 6 and (p (b). In (b), the solid line 
corresponds to the maximum likelyhood, the dashed lines delin- 
eates the one sigma confidence region, and the dotted line, the 
two sigma region. The small circle, at longitude 30° and aper- 
ture 30°, sho ws the approximate lo cation of the hottest region of 
HD 189733b (Kn utson etafll2007h . 



5. Conclusion 

In this study, we provide a simple and very general analytical 
expression of orbital migration induced by anisotropic evapora- 
tion. The amplitude of this effect is not negligible. Indeed, the 
order of magnitude of the migration deduced from the expected 
ejection speed V es and the typical orbital velocity vo of a 3 -day 
orbit is the same as the difference between the observed median 
period of hot Jupiters and hot Neptunes. Moreover, the agree- 
ment is greatly enhanced if ones assumes that the evaporation 
emanates from above the hot test region of the pl anet surface as 
measured on HD 189733b dKnutson et aDl2007l). or predicted 
by atmospheric circulation models ( Showman & Guillot 2002; 
Koskinen et al. 2007; Showman & Polvanil l201ll) . Further stud- 
ies on the dynamics of the exobase may provide stronger con- 
straints on the geometry of evaporation. The analysis presented 
here would then be easily applied to get precise evolution of 
evaporating planets. 
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In order to simplify the study, and to derive analytical ex- 
pressions, we approximate the CDFs of the two populations by 
log-normal distributions 

PDF(lnPi) = Qn Pd = — = expf^^^] , (12) 

1 

where i e {Jup, Nep}, and Pi is expressed in days. N^ (Ti is the 
normal distribution. The best fits give (yUj up = 1.19, <Tj up = 0.39) 
for Jupiter-like planets, and (//Nep = 1.62, <XN ep = 0.57) for 
Neptune-like planets. Assuming that the migration represented 
by A YnP = In 

^Nep — Ift^Jup is independent of Pj U p and PNep> 
the law of total probability implies that A\ nP follows a normal 
distribution with mean //a = //Nep - y^Jup and standard deviation 
°~a = ^"Nep ~ °"jup- Then, a change of variable from A\ nP to r 
Eq. (ITTI) gives 

PDF(r) = - N^ A [-31n(l +r\nq)] . (13) 

1 + r In q 

The PDF (TT31) is plotted in figure Oa) for q = 0.03, which is 
the ratio of the geometric mean of the masses of the two planet 
populations. Using a typical ejection speed V es = 15 km.s -1 and 
the typical initial orbital velocity vo = 145 km.s -1 , the PDF of r 
expressed as a function of (#, <p) constrains the geometry of mass 
ejection (see Fig. Ob)). 

In particular, if one considers an aperture <p lower than 
30°, one obtains the highest probable displacement in the 
eastward direction at about 25° of longitude (solid line in 
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